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1. Introduction 



A staple component of any electrodynamics or quantum mechanics course is the electric po- 
larizability. Neutral material immersed in a weak external field polarizes, internally setting up an 
electric dipole moment, aligned so as to minimize the energy. At the atomic level, the electron 
clouds are distorted creating these microscopic dipole moments. The same process occurs at the 
hadronic level but the polarization effects are now constrained by the strong force. Polarizabili- 
ties of these bound QCD states can be viewed as a distortion of the charged pion cloud of a given 
hadron. 

One can use lattice QCD to non-perturbatively compute the quark and gluon interactions in 
the presence of background electric (or magnetic) fields. For sufficiently weak background fields, 
the low energy properties of the hadrons can be rigorously computed using effective field theory. 
With this treatment, a picture of hadrons emerges from chiral dynamics: that of a hadronic core 
surrounded by a pseudoscalar meson cloud. As some pseudoscalar mesons are charged, polariz- 
abilities of hadrons encode the stiffness of the charged meson cloud (as well as that of the core). 
The form of pseudoscalar meson polarizabilities is consequently strongly constrained by chiral dy- 
namics [||, ||, |J]. However, beyond the leading order, the results depend upon essentially unknown 
low-energy constants, which must currently be estimated in a model-dependent fashion. In the 
case of the charged pion, the experimental measurement of the polarizability has proven difficult, 
both in the original measurement [^] as well as the most recent published result [§]. Currently, 
there is a 2-3 sigma discrepancy between the two-loop %YY prediction and the measured charged 
pion polarizability. New results with higher statistics and the first kaon results are anticipated from 
COMPASS at CERN %. 

A lattice QCD calculation of these quantities is thus very timely. This will allow us to both: 

• make predictions/compare with new experimental results 

• explore stringent predictions from chiral dynamics. For example, the pion polarizabilites 
(from %YY) 

«£ = -^(L 9 + L 10 ), 

o _ a f . s . 
aE ~ 3m n (4xf n y (L1) 



2. Background Field Implementation 

To produce a constant electric field, $ = Si, we use the Euclidean space vector potential, 

A, x (x) = (0,0,-^,0), (2.1) 

where $ is a real-valued parameter. The analytic continuation $ — » — iSm produces a real-valued 
electric field in Minkowski space. This analytic continuation can not generally be performed with 
numerical data due to non-perturbative effects, such as Schwinger pair-creation, which is absent 
in Euclidean space. However, we are only interested in quantities which are perturbative in the 
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.42 



<J>1 = — <J> 2 +2nn, neZ 
(a) (b) 

Figure 1: For a background gauge-potential implemented on a periodic torus (a), the gauge-flux through 
area Al and A2 (b) must be equal and opposite in order to have well defined, single valued action. This 
directly leads to the quantization condition, Eq. (2.3). 



field-strength, for which the naive continuation produces the correct Minkowski-space phsyics, see 
Ref. ^ for explicit details. 

For our initial studies, we implement a quenched electric field (the valence quarks are given 
an electric charge while the sea quarks remain neutral) by post-multiplying existing SU (3) gauge 
links with the Abelian links, 

Ufx (x) -> Up (x) U { / ] (x) , uj? (x) = exp { iQA^ (x) } . (2.2) 

On a torus, construction of a constant field requires quantization to maintain a single-valued ac- 
tion []9[ [H], pH |. The argument is as follows. An action with periodic boundary conditions is defined 
on a torus, which is a closed surface, Fig. [j](a). Therefore, to maintain a single valued action, the 
net flux through the X3-X4 plane, Fig. [[](&) must be zero, modulo 2k, i.e. <J> = QS'fiL = 2%n, with 
nSZ, immediately yielding the 't Hooft quantization condition 

/--^. (23) 

Use of the down quark charge q c j = —\/3e, will lead to the up quark, with charge q u = 2/3e 
automatically satisfying Eq. (|2J|). 



The implementation of a constant gauge field must be slightly modified for a discrete torus Q12| , 



13, 14]. See Refs. ^ [l5j, 16] for our explicit implementation as well as a quantitative study the 
effects on the pion spectra from using non-quantized field values. 



3. Lattice Calculation 

For our calculation, we have used the 2 + 1 flavor 20 3 x 128 dynamical anisotropic gauge en- 
sembles with a renormalized anisotropy of E, = a s /a t = 3.5 and a pion mass of m n ~ 390 MeV JT7| ] 
generated by the Spectrum Collaboration. We used 200 configurations, spaced by either 20 or 40 
monte-carlo time units. On each configuration, for each background field strength, we computed 
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an average of 10 up, down and strange propagators with random spatial sources using the recently 



developed eigCG deflation algorithm Q18|]. For this work, we used the background field strengths 



given by Eq. ( |2.3| ), for n = 0, ±1, ±2, ±3, ±4. On our lattices, this yields the expansion parameter 



governing the deformation of a hadron's pion cloud [|SJ] 

=0.18n 2 , (3.1) 

m\) 

leading us to expect the need of including terms beyond quadratic order in $ to describe the modi- 
fied hadron energies. 

4. Charged Particle Correlation Functions in Background $ -fields 

For neutral particles, it is straightforward to calculate the electric polarizability using standard 
spectroscopy Jl9[]. The lattice two-point functions have their standard forms with the overlap fac- 
tors and energies now depending upon the external field strength, £ . For charged particles, the 
correlation functions are more complicated. Consider a relativistic scalar particle of charge Q. As 
the particle is composite, there are both Born and non-Born terms in the action. The non-Born 
terms account for non-minimal couplings of the field to the particle, such as the polarizabilities. 
These couplings can be summed, as in the case of neutral particles, into the energy of the particle, 
E(S'). 1 Additionally, to describe the two-point correlation functions, one must also sum the Born 
couplings, arriving at ^ 

G{t,g) =Z{g)D{t,E(g),g) +Z'{g)D{t 1 E'{£), + -.., (4.1) 

where 



D(t,E(£),£)= f dsj- — , %-r - . exp 
v v h ' Jo V 27rsinh(£Ms) F 



^— coth(G^) y — L - 



(4.2) 



and 

E{g) =M+-Ana E S' 2 - -^%fa EEE S A + ... (4.3) 



For sufficiently weak fields, Eq. (J4.2 ) reduces to that in Ref. 



D[t,E(*),*) ^ exp l-E{f)t - ^y 4 (^-{E{S)tf + \{E{S)tf + \{E{m) } + 



(4.4) 



5. Results 



In this section, we first present the results of our charged pseudo-scalar mass determinations 
from the long-time behavior of the two-point correlation functions calculated in the presence of 
background electric fields. We then present the resulting determination of the 71* and electric 



polarizabilites, determined from a fit of our extracted masses to Eq. (4.3). The anisotropy of the 



'The quantity E(S) no longer has the interpretation of the energy of the hadron, but remains given by Eq. (|4.3|). 
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Figure 2: Example effective mass plots of charged pion correlation functions, a t M e ft = \n(G(t) / G(t + 1)). 



ensembles used provided a sufficiently fine time resolution as to allow for a fit to the correlation 
functions using a two-state fit, which in turn reduced the systematic uncertainty arising from the 
choice of fitting window. In Fig. ||, we display sample fits to the charged pion and kaon correlation 
functions. The blue data points are the resulting effective masses, 

a t M e& (t)=hx(-^-) . (5.1) 



The yellow band (spanning t = 8 — 29) is the 68% confidence interval using the two-state fit with 
Eqs. (4.1) and (4.2). The gray band is the reconstructed (effective mass of the) correlation func- 
tion using only the ground state correlator and the horizontal salmon colored band is the resulting 
ground state energy, E(£) which enters Eq. (4.2). 

In the upper section of Table |IJ we collect the results of the ground state energies determined 
from this analysis. To enforce the invariance under parity transformations, for each n (from £ = 
Inn/qdfiL) we have taken the average energy from n = ±1, ±2 etc. In the bottom half of Table [j], 
we present the resulting polarizabilites extracted from fitting the pseudo-scalar energies to Eq. (|4~3|). 
In Fig. ||, we present the resulting energy dependence as a function of the background field, and the 
resulting 68% confidence band. The first fit (I) utilizes all five valued of £ while the second fit (II) 
results from a fit to the smallest four background field strengths. 



6. Ongoing and Future Work 



The pion polarizabilites are expected to have significant pion mass dependence, Eq. (1.1). Ad- 
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Table 1: Summary of fit results for charged meson two-point functions for 8 < t < 29. 





n 


a,E{£) 


l-P 


n 




l-P 


K+ 





0.0691(4) 


0.66 


K+ 


0.0969(3) 


0.70 




1 


0.0702(6) 


0.46 


1 


0.0979(4) 


0.77 




2 


0.0718(8) 


0.61 


2 


0.0982(7) 


0.78 




3 


0.0733(16) 


0.93 


3 


0.0958(10) 


0.98 




4 


0.0497(129) 


0.97 


4 


0.0927(23) 


0.97 



71+ 


a t M 


™latt 


Tylatt 

W EEE 


l-P 


I 
II 


0.0692(2) 
0.0692(2) 


18(4)(6) 
16(3) (3) 


24(10) 
17(10) 


0.30 
0.64 


K+ 


a t M 


™latt 


Tylatt 

W EEE 


l-P 


I 
II 


0.0971(2) 
0.0969(2) 


8(3)(1) 
16(4)(3) 


17(5) 
40(9) 


0.03 
0.23 



ditionally, because of this chiral sensitivity, the polarizabilites are also expected to be particularly 
sensitive to the long-range finite volume effects [21]. For these reasons, it is both interesting and 
crucial to perform these calculations with larger spatial volumes and lighter pion masses. This will 
allow for stringent tests of the predictions of #PT as well as first-principles predictions of these 
staple hadronic quantities. 

We are performing similar calculations for the spin- 1/2 baryon polarizabilities 



Acknowledgments 



This work was supported in part by the US Department of Energy. The numerical calculations 



were performed with the Chroma Software Suite [ |23| ] and were carried out on the Jefferson Lab 
High Performance Computers. 



tt + :I 



/r + :II 



0.075 
0.070 
UJ 0.065 
53 0.060 
0.055 
0.050 

0.000 0.005 0.010 0.015 0.020 0.025 0.030 





000 0.005 0.010 0.015 0.020 0.025 0.030 



a t a s e £ a t a s e£> 
Figure 3: The fits to the energy are to all values of n (I) and up to n — ±3 (II). 1 — P is the quality of fit. 
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